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Autophagy is a catabolic pathway typically induced by nutrient starvation to recycle amino acids, but can also func- 
tion in removing damaged organelles. In addition, this pathway plays a key role In eukaryotic development. To date, not 
much Is known about the role of autophagy in apicomplexan parasites and more specifically in the human malaria para- 
site Plasmodium falciparum. Comparative genomic analysis has uncovered some, but not all, orthologs of autophagy- 
related (ATG) genes in the malaria parasite genome. Here, using a genome-wide In silico analysis, we confirmed that ATG 
genes whose products are required for vesicle expansion and completion are present, while genes involved In Induction 
of autophagy and cargo packaging are mostly absent. We subsequently focused on the molecular and cellular function 
of P. falciparum ATG8 (PfATGS), an autophagosome membrane marker and key component of the autophagy pathway, 
throughout the parasite asexual and sexual erythrocytic stages. In this context, we showed that PfATGS has a distinct 
and atypical role in parasite development. PfATGS localized in the apicoplast and In vesicles throughout the cytosol dur- 
ing parasite development. Immunofluorescence assays of PfATGS In apicoplast-mlnus parasites suggest that PfATGS is 
involved In apicoplast biogenesis. Furthermore, treatment of parasite cultures with bafllomycin A^ and chloroqulne, both 
lysosomotropic agents that Inhibit autophagosome and lysosome fusion, resulted In dramatic morphological changes 
of the apicoplast, and parasite death. Furthermore, deep proteomic analysis of components associated with PfATGS Indi- 
cated that It may possibly be Involved In ribophagy and piecemeal microautophagy of the nucleus. Collectively, our 
data revealed the importance and specificity of the autophagy pathway in the malaria parasite and offer potential novel 
therapeutic strategies. 



Introduction 

Malaria continues to be a major global health problem with 
up to a million deaths per year.' The human malaria parasite 
has a complex life cycle consisting of 3 major stages in the mos- 
quito, the human liver, and the human blood. In the blood stage, 
merozoites invade erythrocytes and within 48 h the parasite pro- 
gresses through ring, trophozoite, and schizont stages to produce 
up to 32 new daughter cells that egress and reinvade new red 
blood cells. This asexual cycle continues until the human host 
is stressed. Parasites then undergo gametocytogenesis, which 
can require up to 2 wk for full differentiation. When another 
female mosquito takes a blood meal, mature male and female 
gametocytes are taken up and undergo sexual reproduction in 
the mosquito midgut to complete the cycle. Throughout its life 
cycle, the parasite undergoes cell differentiation and metamor- 
phosis multiple times, requiring degradation of unnecessary 
stage-specific proteins and organelles within host cells. In the 
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face of increased resistance of P. falciparum to current antima- 
larials, understanding the mechanisms that control the parasite 
degradation pathways may open new doors for the development 
of novel antimalarial drugs. 

In eukaryotic cells, the autophagosome-lysosome pathway is 
involved in degradation and recycling of proteins and organelles. 
Similar to the ubiquitin-proteasome system that degrades short- 
lived and misfolded proteins,^ the autophagy pathway consists 
of a number of sequential steps: selection and tagging of cargo, 
recognition and transport to proteolytic machinery, degradation 
of cargo, and recycling of small polypeptides or amino acids. 
Proteins and organelles targeted for degradation are enveloped 
by a double membrane to form an autophagosome that fuses 
with a lysosome to degrade the cargo. There are 3 main mem- 
brane-mediated processes conserved through most eukaryotes: 
1) macroautophagy (referred to as autophagy hereafter) removes 
unnecessary proteins and damaged organelles, 2) microautoph- 
agy occurs when the lysosome directly envelops cytoplasmic 
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material, and 3) chaperone-mediated autophagy (CMA) relies on 
a chaperone protein to translocate individual unfolded proteins 
into the lysosome. Once thought to be a nonspecific molecular 
process that recycles proteins for survival under nutrient-lim- 
ited conditions, autophagy is involved in aging, degradation of 
pathogens, cell differentiation, and removal of old or damaged 
organelles.^ 

More than 30 autophagy-related {ATG) genes were first dis- 
covered in yeast,'' and were thought to be generally conserved in 
eukaryotes.' However, preliminary data show that the canonical 
machinery of autophagy in parasitic protists has become diverse 
as a possible consequence of their adaptive life cycles.'' ylTTJ genes 
have been modified, lost or have expanded, but, overall, auto- 
phagy seems to play a critical role in parasite development.^ For 
example, in Toxoplasma gondii, an apicomplexan parasite related 
to Plasmodium, orthologs of the ATG12-ATG5 conjugation sys- 
tem have not yet been identified. However, the ATG8 conjuga- 
tion system has been detected and ATGl and TOR orthologs 
have been found to mediate autophagosome formation upon 
starvation.** In addition, it was recently found that in Toxoplasma, 
TgATG3, an E2 conjugating enzyme for the LC3 lipidation pro- 
cess, is essential for regulating mitochondrial homeostasis during 
tachyzoite development.' In the murine rodent malaria parasite, 
P. berghei, autophagy is required to remove unnecessary organ- 
elles during metamorphosis of the sporozoite into merozoites in 
hepatocytes.'" 

Despite the publication of these recent studies, not much is 
known about the role of autophagy in the human malaria para- 
site, P. falciparum. Genome-wide analyses using BLAST searches 
for either homologs or putative orthologs of known autophagy 
components have uncovered a rudimentary set of ATG genes. 
Recently, PfATGS, a marker of the phagophore and auto- 
phagosome membranes, was shown to localize in the apicoplast 
membrane at the late schizont and merozoite stages." In the rodent 
malaria, P. berghei, the function of PbATGS seems restricted to 
apicoplast biogenesis in the liver stage. The apicoplast is a plas- 
tid of secondary endosymbiotic origin that is essential to parasite 
survival.'^ The role and function of PfATGS in the apicoplast 
remains ambiguous, but it has been implied that PfATGS may be 
part of endoplasmic reticulum (ER) -related organelle biogenesis, 
which includes apicoplast membrane biogenesis.''' In this study, 
we used a combination of genomic, molecular, cellular, and pro- 
teomic approaches to further investigate the role oi ATG genes 
throughout the P. falciparum asexual and sexual blood stages. 
We demonstrated that the human malaria parasite has modified 
the canonical autophagy pathway to a rudimentary set oi ATG 
genes that participate not only in apicoplast biogenesis, but in 
additional functions of cellular development. 

Results 

Identification of autophagy genes in P. falciparum 

Previous in silico analysis of P. falciparum's genome using 
BLAST retrieved few homologs of ATG genes. To uncover a 
complete list oiATG genes in P. falciparum, we conducted a hid- 
den Markov model (HMM) search for ATG protein domains. 



The HMM methodology has a high sensitivity for the detection 
of less conserved regions, and is therefore better applicable to api- 
complexan parasites which seem to have modified the canonical 
autophagy pathway.'^ 

Of the upstream autophagy pathway, none of the members 
of the Atgl kinase complex functional unit were detected, even 
though putative Atgl and Atgl7 orthologs have been reported 
previously (Table SI)." Likewise, the transmembrane protein 
Atg9 was not identified. However, we detected potential candi- 
dates for PfATG6 and PfATGl4, members of the class III Ptdlns 
3-kinase complex. In addition, a possible PfATG2 candidate with 
homology in the C-terminal region was identified.'^ Altogether, 
this in silico study suggests that the human malaria parasite lacks 
proteins critical for detection of nutrient limitation and induc- 
tion of the canonical autophagy pathway, unlike higher eukary- 
otes and the apicomplexan T. gondii. 

Of the downstream pathway, we identified most genes in 
the Atgl 2 conjugation system with the exception of the E2-like 
enzyme AtglO (Table SI). In addition, all members of the AtgS 
conjugation system (Atg7, Atg3, and Atg4) were detected. 
Finally, a homolog of the lipase Atgl5, responsible for breakdown 
of intravacuolar membranes, was identified."* As a whole, our 
comparative genomics analysis has expanded the ATG gene list 
found in P. falciparum and indicates that a rudimentary auto- 
phagy pathway is present and likely functional. 

PfATGS and GFP-PfATGS expression and localization dur- 
ing the asexual erythrocytic stage 

To investigate the role of autophagy in the human malaria 
parasite at the molecular and cellular levels, we first explored the 
role of ATGS during the erythrocytic cycle. ATGS is the stan- 
dard autophagosome marker and is essential to autophagosome 
formation. It remains present on the outer autophagosome mem- 
brane until, or just before, fusion with the lysosome, and a popu- 
lation of ATGS also remains bound to the inner autophagosome 
membrane." We generated a custom PfATGS-specific antibody 
and analyzed PfATGS protein expression in ring, trophozoite, 
and schizont stage parasites. Protein concentrations of parasite 
lysates from each stage were determined by Bradford assay^° and 
equal amounts of protein were resolved by SDS-PAGE and pon- 
ceau staining (Fig. SIA). PfATGS expression was detected at all 
stages of the asexual cycle, with a decreased abundance at the 
trophozoite stage (Fig. lA; Fig. SIA). Similar to previous reports 
that suggest PfATGS is primarily present in a phosphatidyletha- 
nolamine (PE) -conjugated form," we observed PfATGS as a 
single band (Fig. lA; Fig. SIA). In most other eukaryotes, ATGS 
requires cleavage by ATG4 to expose a glycine for conjugation to 
PE. In contrast, the terminal residue of PfATGS is a glycine and 
we therefore speculate that it is predominately lipidated during 
the asexual cycle. 

Subcellular fractions of mixed-stage parasites were immunob- 
lotted for PfATGS. Examination of cytoplasmic, nuclear, and 
membrane fractions in the 3D7 cell line indicated that PfATGS 
was predominantly incorporated into membranes (Fig. IB; 
Fig. SIB). To validate our fractionation, we used the AGP 
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GFP cell line, which expresses GFP in the cytosol, as a blotting 
control for the cytosolic fraction (GFP) and the nuclear fraction 
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Figure 1. Expression of PfATGS throughout the erythrocytic cycle. (A) Immunoblot analysis of 
PfATGS expression in tightly synchronized parasites harvested at the ring, trophozoite, and schiz- 
ont stage. Protein concentration was determined by Bradford analysis and 30 |j,g of parasite lysate 
was loaded in each lane. Using a custom-generated anti-PfATGS antibody, it was observed that 
PfATGS protein is present throughout the asexual cycle, albeit at lower abundance in the tropho- 
zoite stage. (B) Immunoblot analysis showing the expression of endogenous PfATGS in cytosolic, 
nuclear, and membrane fractions of parasite lysates. Mixed-stage 3D7 parasite cultures were used 
for subcellular fractionation. Protein concentration was determined by Bradford assay and 30 [j,g 
of protein was loaded for the cytosolic and nuclear fractions, while 10 [lq of protein was loaded 
for the membrane fraction. PfATGS is detected in all fractions, but is predominantly membrane- 
bound. (C) Immunoblot analysis of the GFP-PfATGS transiently-transfected P. falciparum cell line. 
IVIixed-stage 3D7 and GFP-PfATGS parasites were lysed and 30 p,g of lysate was loaded into each 
lane. Endogenous PfATGS was detected using anti-PfATGS at 14 kDa, while the fusion protein GFP- 
PfATGS was detected by both anti-GFP and anti-PfATGS antibodies at the expected molecular 
weight of 41 kDa. For all panels, relative amounts of protein loaded were assessed by Ponceau- 
staining the full-length PVDF membrane prior to blotting (Fig. SI). The sizes of molecular mass 
markers are indicated in kDa. 



(histone H3) (Fig SIC). To further define the subcellular local- 
ization of PfATGS, synchronized 3D7 parasite cultures were 
immunolabeled with anti-PfATGS at all stages of the asexual 
cell cycle. PfATGS was found on punctate structures of 200- 
to 350-nm diameter throughout the asexual cycle (Fig. S2). 
These structures are smaller than the classical autophagosome 
of 300- to 900-nm diameter and larger than the typical 150- 
nm diameter cytoplasm-to-vacuole targeting vesicles seen in 
yeast. During the ring stage, small punctate structures were 
observed, which merged into a circular structure during the tro- 
phozoite stage with additional puncta observed in the cytosol of 
P. falciparum and host red blood cell (RBC; Fig. S2). This circu- 
lar structure branched out as the parasite underwent schizogony, 
and finally segmented into daughter cells (Fig. S2). To comple- 
ment our results, we generated a GFP-PfATGS transiently trans- 
fected cell line that expresses the fusion protein at the expected 



molecular weight of 47 kDa (Fig. IC; 
Fig. SID). No apparent difference in GFP- 
PfATGS parasite growth rate or propaga- 
tion was observed, although blood smears 
of parasites did have an enlarged digestive 
phagolysosome (food vacuole) pheno- 
type at the trophozoite stage (Fig. S3B). 
The GFP-PfATGS-transfected cell line 
displayed a similar pattern as compared 
with anti-PfATGS labeling with a diffuse 
GFP signal throughout the cytoplasm 
(Fig. S4). Immunoblot analysis showed 
that while GFP-PfATGS was also mem- 
brane-bound, GFP-PfATGS was more 
abundant than endogenous PfATGS in 
cytosolic fractions (Fig. ID; Fig SIE). 
Interestingly, at the trophozoite stage the 
localization of GFP-PfATGS and PfATGS 
exhibited less overlap than at other stages. 
Metabolic activity is high during the tro- 
phozoite stage and as vesicle transport 
increases, the exogenous PfATGS and 
GFP-PfATGS may be recycled onto newly 
synthesized vesicles. PfATGS was recently 
reported to be solely localized in the 
apicoplast during the merozoite stage." 
However, we observed PfATGS-decorated 
vesicles, presumably autophagosomes, at 
all stages of the asexual blood cycle, as 
recently shown by Tomlins et al.^' 

PfATGS vesicles and localization 
with the apicoplast 

The transgenic DIO strain with the 
apicoplast acyl carrier protein (ACP) fused 
to GFP was used to confirm the possible 
colocalization of PfATGS and the apico- 
plast.^" PfATGS was indeed detected close 
to the apicoplast, suggesting a possible 
biological role of this protein in the bio- 
genesis of the organelle (Fig. 2) . However, 
additional small punctate structures that do not completely over- 
lap with the apicoplast were observed throughout the asexual 
cycle. Quantification of vesicles revealed 1 to 2 puncta at the ring 
stage, while 10 to 15 puncta, and 15 or more were observed at the 
trophozoite and the schizont stages, respectively. Moreover, the 
branching of PfATGS during the trophozoite and schizont stages 
covered a larger area than ACP-GFP and appeared bulbous com- 
pared with the smooth branching of the apicoplast. In this para- 
site line, vesicles decorated with PfATGS were also observed in 
the host erythrocyte from the late ring to late trophozoite stages. 

Recently, it was reported that apicoplast-minus P. falciparum 
strains could be established by treatment of the parasites with 
antibiotics that inhibit apicoplast gene transcription and transla- 
tion, followed by rescue of the parasites by isopentenyl pyrophos- 
phate (IPP) supplementation.^^ In antibiotic-treated ACP-GFP 
parasites rescued with IPP, apicoplast proteins were localized in 
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Figure 2. PfATGS immunofluorescence 
staining throughout the erythrocytic 
cycle. Tightly synchronized P. falciparum 
D10-ACP-GFP parasites were immuno- 
labeled with anti-PfATG8. D10 parasites 
express ACP-GFP, an apicoplast signal 
and transit peptide fused to GFP that 
localizes to the apicoplast. Localization 
of ACP-GFP and PfATGS can be observed, 
although PfATGS appears to surround the 
apicoplast, indicating a possible associa- 
tion with the outer membrane. Additional 
PfATGS vesicles can be observed in the 
cytosol and host red blood cell from the 
late ring to late trophozoite stage, when 
parasites are metabolically active. DAPI is 
used as a DNA marker in blue. Scale bar: 
5 |jLm. 



numerous small foci spread through- 
out the cell with additional diffuse 
GFP staining, while branching of the 
apicoplast at the schizont stage was 
completely abolished. Labeling of api- 
coplast-minus parasites with PfATGS 
showed that PfATGS vesicles overlap 
with some, but not all, ACP-GFP vesi- 
cles (Fig. 3). Quantification of vesicles 
revealed that two-thirds of PfATGS 
vesicles did not overlap with the ACP- 
GFP foci. 

PfATGS localization during 
gametogenesis 

Next, we investigated the distri- 
bution of PfATGS during the sexual 
phase of the erythrocytic stage. In con- 
trast to the asexual cycle which is com- 
pleted within 48 h, differentiation into 
gametocytes is a relatively slow process 
requiring 8-10 d.^** In the NF54 strain, 
that readily produces gametocytes, 
the apicoplast ACP-GFP cell line and 
our newly generated GFP-PfATGS 
parasite line, gametocytes of all stages 
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Figure 3. PfATGS immunofluorescent staining in 
apicoplast-minus malaria parasites. The integrity of 
the apicoplast in these parasites was lost using treat- 
ment with the antibiotic chloramphenicol. GFP-ACP 
parasite cultures were rescued by supplementing 
media with isopentenyl pyrophosphate (IPP), the 
essential product from the non-mevalonate isopren- 
oid precursor pathway in the apicoplast. The large 
apicoplast PfATGS-labeled branching, as observed 
in untreated parasites, is disrupted into multiple foci 
in apicoplast-minus parasites. These foci show colo- 
calization of GFP-ACP and PfATGS vesicles, although 
additional PfATGS vesicles can be observed in the 
cytosol. DAPI is used as a DNA marker in blue. Scale 
bar: 2.5 /jum. 
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Figure 4. PfATGS immunofluorescent staining in gametocytes. (A) Validation of PfATGS localization 
using the ACP-GFP cell line. A clear localization of GFP-PfATGS and the apicoplast can be observed. 

(B) Expression of PfATGS in GFP-PfATGS gametocytes. An abundance of PfATGS-labeled autophagosomes 
is observed during gametocytogenesis. In stage ll-IV, large vesicles are observed at the apical poles. 

(C) Differential PfATGS pattern in male and female gametocytes. Female NF54 gametocytes at stage IV 
have rounded poles and PfATGS vesicles are more uniform in size and distribution in the cytosol, while 
male gametocytes have larger PfATGS vesicles at the apical poles (white arrowheads). DAPI is used as a 
DNA marker in blue. Scale bar: 5 ixm. 



showed numerous PfATGS-associated vesicles throughout the 
cytosol. In the ACP-GFP cell line, it was observed that PfATGS 
localized in the apicoplast adjacent to the nucleus (Fig. 4A). 



During stage II-IV, PfATGS was also 
present at the apical poles (Fig. 4B). 
Interestingly, the apical pole stain- 
ing differed depending on the sex of 
gametocytes. Stage IV male gameto- 
cytes displayed a large vesicle at both 
poles, whereas female gametocytes 
had vesicles more evenly dispersed 
outside the nucleus (Fig. 4C). Once 
gametocytes are in the mosquito gut, 
females will produce a single macro- 
gamete whereas males will produce 
8 microgametes. During microgame- 
togenesis, parasites exflagellate and 
develop axonemes, but the apicoplast 
and mitochondrion remain in the 
original gametocyte cytoplasm.-^ A 
temporal and sex-dependent regu- 
lation of the autophagy pathway 
in male and female gametocytes is 
therefore not surprising. 

Absence of induction of the auto- 
phagy pathway by starvation in 
P. falciparum 

Autophagy is typically induced by 
nutrient starvation.* However in P. 
falciparum, the ATG genes involved 
in autophagy induction were not 
detected in our in silico genome-wide 
analysis (Table SI). To validate our 
in silico data and test whether auto- 
phagy could be induced by nutrient 
limitation in P. falciparum, parasite 
cultures were grown in both serum- 
free medium and minimal medium 
devoid of the essential amino acid 
L-glutamine and serum, with mini- 
mal glucose. Ring stage parasites 
grown in medium without human 
serum or albumax for 24 h were 
arrested at the trophozoite stage, 
and normal growth was restored 
upon addition of complete medium. 
A similar state of hibernation has 
been recently observed in isoleucine- 
starved parasites.^'' Parasites that 
were starved for 72 h in minimal 
medium and then given complete 
medium produced 1% parasitemia 
of gametocytes when compared with 
control parasites. However, despite 
multiple attempts, incubation in 
minimal medium for 6, 24, or 72 h 
did not increase the amount of PfATGS autophagosomes (data 
not shown). These data indicate that unlike other eukaryotic 
cells, autophagy in the asexual blood stage of P. falciparum 
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Figure 5. Effect of antimalarial treatment on PfATGS expression and localization. (A) PfATGS immunofluorescent staining in 3D7 parasites treated with 
lysosomotropic agents and artemisinin. After 3 h incubation with 75 nM of bafilomycin A,, branching of the apicoplast is disrupted and multiple PfATGS 
autophagosomes can be observed. Artemisinin treatment at a high concentration (25 (jlIVI) results in cell cycle progression arrest, and chloroquine treat- 
ment, at a concentration of 100 nM (IC^^), causes disruption of the apicoplast integrity. DAPI is used as a DNA marker in blue. Scale bar: 2.5 ^im. (B) Graph 
of PfATGS vesicles on trophozoite-stage parasites with or without autophagy inhibitor incubation (baf, 75 nM bafilomycin A,; art, 25 jj,M artemisinin; cq, 
100 nM chloroquine). (C) Western blot analysis of PfATGS expression levels upon bafilomycin A,, chloroquine and artemisinin treatment. 



cannot be activated under mild starvation conditions, and that 
parasites starved for long periods of time preferably undergo 
gametocy togenesis . 

Downstream inhibitory effects on P. falciparum's auto- 
phagy pathway 

Traditional methods to induce and inhibit autophagy do not 
seem to be effective in P. falciparum. The lack of classical induc- 
tion by starvation or drug treatments may be due to the parasite's 



divergent upstream signaling kinases. Indeed, treatment of para- 
site cultures with the TOR complex inhibitor rapamycin did not 
induce autophagy (data not shown), most likely because the para- 
site lacks a homolog for TOR. Similarly, PfATGS localization was 
previously reported to remain unchanged upon treatment of para- 
sites with wortmannin, a PtdIns3K and class I PI3K inhibitor." 
Both rapamycin and wortmannin work by targeting upstream 
proteins of the autophagy pathway. Bafilomycin A^, on the other 
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hand, is a V-type H*-ATPase inhibitor that inhibits the down- 
stream fusion events between autophagosomes and lysosomes, 
and has effects on transport within the endocytic pathway.^^ 
Here, we observed that bafilomycin is a potent inhibitor of the 
chloroquine-sensitive 3D7 strain, and the chloroquine-resistant 
strains W2 and Dd2, as well as our GFP-PfATG8 cell line, with 
IC^j, values ranging between 21 to 26 nM (Table S2). The mor- 
phological effects of bafilomycin A^ were examined by treatment 
of synchronized 3D7 and GFP-PfATG8 parasite cultures at the 
ring, trophozoite, and schizont stages for 6 and 12 h. A rapid inhi- 
bition of parasite growth was observed, in particular at the mature 
stages (Fig. S3). The potent inhibitory effect of bafilomycin A^ 
on parasite growth could be reversed by washing bafilomycin A^ 
out of the culture medium within 3 h (data not shown). After 3 
h of bafilomycin Aj incubation at 75 nM, early trophozoite para- 
sites decreased in size with a slight decrease of PfATGS -labeled 
vesicles, from an average of 11 punctate structures to an average of 
8.5 puncta {P = 0.005) (Fig. 5A and B). This overall decrease of 
PfATG8 staining was accompanied by dissociation from the api- 
coplast and an increased staining level in the cytosol. Immunoblot 
analysis of PfATGS expression after bafilomycin A^ treatment 
did not show a substantial difference when compared with the 
control (Fig. 5C; Fig. SIF). Moreover, immunoblot analysis of 
cytoplasmic and membrane fractions did not indicate changes 
in subcellular PfATG8 expression (data not shown). Similar to 
what we observed at the trophozoite stage, schizont-stage parasites 
incubated with bafilomycin A^ displayed a strong and significant 
dissociation between PfATGS and the apicoplast, and an increase 
in PfATGS vesicle staining from an average of IS to an average of 
30 puncta {P = 0.0003) (Fig. 5A and B). 

In eukaryotic cells, autophagy plays a role in a particular form 
of cell death. P. falciparum was previously reported to undergo 
vacuolization after 24 h of incubation with antimalarial com- 
pounds."** We examined the effects of artemisinin, currently the 
most effective antimalarial agent,^' and chloroquine on PfATGS 
localization. ACP-GFP parasites treated at the trophozoite stage 
with a high concentration (25 |xM) of artemisinin did not show 
any significant changes in PfATGS staining pattern (Fig. 5A). 
Parasites did exhibit a delay in cell cycle progression and 
remained stunted at the trophozoite stage when artemisinin was 
added (Fig. 5A). On the other hand, treatment with chloroquine, 
which was once the most successful antimalarial drug therapy, 
but is also often used as a tool to study the autophagy pathway in 
higher eukaryotes, showed a strong phenotype. In all eukaryotic 
cells, chloroquine increases the lysosomal pH and prevents the 
fusion of autophagosomes with lysosomes, and therefore blocks 
lysosomal protein degradation. Trophozoite parasites treated 
with 100 nM of chloroquine (ICg^) for 6 h displayed fragmen- 
tation of the apicoplast with no significant change in PfATGS 
vesicles from an average of 15.5 (in the absence of chloroquine) to 
an average of 17.4 puncta (with chloroquine treatment) (Fig. 5B). 
Whether this phenotype is linked to vesicle traffic inhibition or 
an active fragmentation process remains to be determined. 
Identification of proteins associated with PfATGS 
To further confirm the role of the autophagy pathway in para- 
site development and to identify additional proteins involved in 



this process, we performed an immunoprecipitation experiment 
using an anti-GFP antibody on synchronized GFP-PfATGS 
cultures in combination with a shotgun mass spectrometry 
approach. Proteins were identified and quantified using multi- 
dimensional protein identification technology (MudPIT) cou- 
pled with spectral counting.'" Parasites of early trophozoite and 
late schizont stages were either nontreated or were treated with 
bafdomycin A^ for 3 h before harvesting to potentially enrich for 
PfATGS-associated proteins. Proteins were considered positive 
hits if they passed 2 sets of filtration criteria. First, proteins had 
to be detected in at least 2 of the 3 technical replicates of at least 
one sample type, and had to be either not detected in our 2 inde- 
pendent negative controls, or be more than 5 -fold enriched com- 
pared with these negative controls (a precipitation without the 
anti-GFP primary antibody, and a GFP-precipitation from the 
3D7-HT-GFP cell line that stably expresses GFP in all stages of 
the parasite life cycle).'' Second, proteins also had to be detected 
in at least 4 runs (out of a total of 12 runs) and not in controls, 
or if detected in controls with a QSpec z-score'^ above 1.7 in at 
least 1 of the 4 sample types (see Materials and Methods). The 
lack of significant overlap between proteins from our precipita- 
tions and controls validates the strength and specificity of our 
approach. In total, 101 proteins associated with GFP-PfATGS 
vesicles passed our selection criteria (Table S3). The autophagy 
network in human cells undergoing basal autophagy consists 
of 6 ATGS orthologs that interact directly with 67 proteins.'' 
If the autophagy pathway indeed has multiple functions in the 
malaria parasite, this number of proteins associating with ATGS 
vesicles is not surprisingly high. Evidence for enrichment of 
PfATGS-associated proteins by bafilomycin Aj treatment was 
observed in the schizont stage, where we detected 77 proteins in 
at least 2 replicates of the treated sample vs. 46 such proteins in 
the untreated sample. This correlates with an increased amount 
of PfATGS vesicles that were observed when parasites were 
treated with bafilomycin A^. Among the proteins found to be 
associated with PfATGS were PfATG7 (El-like) and 2 undistin- 
guishable isoforms of PfATG3, whereas PfATG6 was only sig- 
nificantly detected in the drug-treated schizont stage. PfATG2 
and PfATGlS were detected in the whole data set, but did not 
pass our selection criteria, whereas PfATG5 was only detected in 
one of our controls. 

Subsequently, a gene ontology (GO) analysis was performed 
to find functional enrichments among the proteins associ- 
ated with PfATGS. A relatively large overlap of significantly 
enriched GO terms was observed between the different samples 
(Fig. S5A). Among all 101 proteins that were found to be asso- 
ciated with PfATGS, significantly enriched GO terms included 
those involved in the digestive phagolysosome, the phagolyso- 
some, protein transport, nuclear import, and metabolic processes 
(Fig. S5B; Table S4). 

In line with the fact that more than half of P. falciparum's pro- 
teome is incompletely characterized, approximately one-third of the 
proteins identified to be associated with PfATGS were hypothetical 
proteins. These hypothetical proteins were BLASTed against the 
yeast genome to find evidence for their function. We found homol- 
ogy of PF3D7_1 339700 to the yeast protein Nvjl that is present 
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at the formation of nucleus-vacuole junctions and is involved in 
piecemeal microautophagy of the nucleus.^'*'^* In addition, we 
detected homology of PF3D7_1459000 to yeast Bre5, which is 
associated with the process of ribophagy.'^'^** Further validation will 
be required to determine if these homologous proteins are indeed 
involved in selective autophagy pathways and whether they have a 
functional role in the human malaria parasite. Collectively, these 
data indicate that a functional PfATGS conjugating system is pres- 
ent at the proteomic level and highlight the multiple facets of this 
pathway in the parasite's infectious cycle progression. 

Discussion 

The data presented in this study provide a basic understand- 
ing of the role of autophagy in the human malaria parasite. 
Our in silico, cellular, and molecular approaches demonstrate 
that the canonical autophagy pathway has become modified in 
P. falciparum. Similar to T. gondii, P. falciparum encodes homo- 
logs for only a subset of yeast ATG proteins. "^'^ Upstream kinases 
involved in detection of nutrient limitation and induction of 
autophagy are absent. P. falciparum is an intracellular parasite 
residing inside the red blood cell where hemoglobin, its primary 
nutritional source, is abundant. The parasite may have lost genes 
required for the nutrient-dependent induction of autophagy, as 
seen in other systems.^' This was further confirmed by our obser- 
vation that autophagy cannot be induced by mild nutrient starva- 
tion, or by specific inhibitors of key upstream kinase proteins in 
the pathway that are involved in regulation in other organisms. 
Our results are in agreement with a recent study that showed 
that P. falciparum enters a state of hibernation in isoleucine-free 
medium.-'' Effector proteins that initiate the autophagy pathway 
in the human malaria parasite remain unknown. 

On the contrary, most downstream proteins involved in auto- 
phagy were detected in P. falciparum and appear to be functional. 
PfATG proteins involved in autophagosome formation, retrieval, 
and vesicle breakdown were clearly identified. From the 12 ATG 
genes uncovered from our HMM search, our proteomic analy- 
sis identified 7 PfATG proteins, although PfATG2, PfATG5, 
PfATG6, and PfATGlS were detected at low abundance levels. 
Proteins of the ATG8 conjugation system (PfATG7, PfATGS, 
and 2 isoforms of PfATG3) were significantly enriched in the 
GFP-ATG8 pull-downs, with the exception of PfATG4. Our gel- 
free mass spectrometry analysis identified the E2-like PfATG3, 
validating previous reports of a structural characterization of 
PfATG8-PfATG3 interaction."" The role of PfATG4 in priming 
PfATGS for PE lipidation remains uncertain. PfATGS was pri- 
marily found conjugated to PE, and proteolytic cleavage by PfAtg4 
may not be required for lipidation since PfATGS's C-terminal 
residue is a glycine. A recent report in T. gondii demonstrated 
that TgATG4 recycles membrane-bound TgATGS,'" but fur- 
ther investigation will be required to determine the exact role of 
PfATG4 in recycling PfATGS from membranes in Plasmodium. 

In addition to the detection of PfATGS localizing with the 
apicoplast"''^ and its possible involvement in apicoplast biogen- 
esis, we report for the first time the presence of PfATGS-coated 
vesicles throughout the entire asexual and sexual erythrocytic 



cycle. Surprisingly, PfATGS vesicles were also observed within 
the host RBC. Mature erythrocytes are denucleated and lack 
proteins involved in vesicle traffic. In order for parasites to trans- 
port parasite proteins to the host RBC membrane, parasites have 
expanded the classical secretory pathway by exporting machinery 
into the RBC that controls vesicle-mediated traffic.''^ Currently 
we can only speculate that PfATGS decorated vesicles in the RBC 
cytosol are part of the P. falciparum exported secretory system. 
Recently, exosome-like vesicles containing nuclear material were 
found to be mediators of intercellular communication between 
parasites.*^ Interestingly, our GO analysis of PfATGS-associated 
proteins revealed enrichment for the nucleus and nuclear enve- 
lope. It is unknown what the contents of PfATGS vesicles are and 
if these structures are transported to Maurer's clefts where cargo 
can be stored for export to the RBC membrane.''* 

Over the years, details of protein trafficking to the apicoplast 
have remained ambiguous. '''•^^'''^'''* Apicoplast proteins derived 
from nuclear-encoded genes utilize a bipartite N-terminal sig- 
nal peptide to enter the secretory pathway and a transit peptide 
for apicoplast import.^^ '" The 2 outer membranes of the apico- 
plast are ER-related.''^ Similarly, the phagophore may also derive 
from the ER membrane.'"*""^ It is therefore possible that these 
2 pathways could interact with each other. In apicoplast-minus 
parasites, PfATGS localized with vesicles that contain apicoplast- 
targeted proteins, which suggests an alternative form of bulk 
protein trafficking to the apicoplast. Further examination into 
PfATGS cargo proteins will be necessary to determine the role 
of PfATGS and the autophagy machinery's involvement in api- 
coplast biogenesis. 

Interestingly, known inhibitors of the autophagy pathway and 
parasite development seem to interfere with apicoplast biogenesis. 
The branch-like structure of the apicoplast observed during the 
schizont stage was destroyed into multiple foci in both apico- 
plast-minus parasites and cultures treated with lysosomotropic 
agents. Both chloroquine and bafilomycin Aj cause vacuolar 
deacidification and inhibit fusion of autophagosomes with lyso- 
somes. Bafilomycin A^ also blocks transport within the endocytic 
pathway," but this is the first report of bafilomycin A^ blocking 
transport to a plastid organelle. Chloroquine is an antimalarial 
that accumulates in the digestive phagolysosome and inhib- 
its the formation of hemozoin, the byproduct of hemoglobin 
catabolism.''^ In human fibroblasts, chloroquine also raises intra- 
lysosomal pH and disrupts the intracellular pathway for newly 
synthesized acid hydrolases and uptake by pinocytosis.''^ Reports 
of cytoplasmic vacuolization, similar to autophagy-associated cell 
death, have been observed in parasites after chloroquine treat- 
ment, but whether or not unicellular protists undergo apoptosis 
remains ambiguous.^'* Our lysosomotropic drug-treated parasites 
displayed an accumulation of PfATGS vesicles and apicoplast 
fragmentation, which could be interpreted as either cytoplasmic 
vacuolization or an inhibition of vesicle fusion. In comparison, 
incubation with the antimalarial artemisinin did not induce 
autophagy-associated cell death or vacuolization. Further investi- 
gations will be needed to determine whether these lysosomotropic 
agents inhibit fusion with the apicoplast or cause apicoplast frag- 
mentation by vacuolization. 
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While the asexual stage of the erythrocytic cycle does not last 
for more than 48 h, the differentiation into gametocytes can take 
up to 2 weeks.''' During this time, the parasite resides within 
the same red blood cell and will need to recycle its proteins in 
order to survive. The abundance and atypical localization of 
PfATGS in gametocytes suggests that autophagy indeed plays an 
important role in protein turnover and parasite survival. It has 
been previously reported that the autophagy pathway degrades 
unnecessary organelles, such as micronemes, during sporozoite 
metamorphosis into merozoites in the liver stage.'" Our data 
show that localization of PfATGS remains on the apicoplast, but 
there is an increase in autophagosomes during gametocytogen- 
esis. Considering the apparent role of PfATGS in gametocyte 
development, the observation that harsh nutrient starvation 
induced gametocytogenesis within 48 h without reinvasion of a 
new erythrocyte may be indicative of a specialized role of auto- 
phagy in cell differentiation. Possibly, a delayed or modified 
form of autophagy induction could trigger differentiation into 
the sexual stage in P. falciparum. Further studies into this process 
are required to determine the developmental cues necessary for 
gametocytogenesis. 

Collectively, this study presents the critical role of a modi- 
fied autophagy pathway in the human malaria parasite. PfATGS 
has evolved to participate in multiple functions including, 
but not limited to, cell differentiation, protein turnover, 
apicoplast biogenesis, and partitioning into daughter cells. 
While macroautophagy is the best-characterized type of 
autophagy, our proteomic approach validates the importance 
of this pathway, but suggests that ribophagy and piecemeal 
microautophagy of the nucleus may also occur in P. falciparum, 
although additional experiments will be necessary to validate 
these findings. A detailed understanding of the molecular mech- 
anisms regulating this pathway throughout the human malaria 
parasite life cycle will not only increase our knowledge of para- 
site biology, but will most likely uncover novel antimalarial 
strategies. 

Materials and Methods 

P. falciparum parasite culture 

The following reagents were obtained through the MR4 as 
part of the BEI Resources Repository, NIAID, NIH: Plasmodium 
falciparum iVfl, MRA-102, deposited by DJ Carucci; Dd2, 
MRA-156, deposited by TE Wellems; DIO ACP (leader) -GFP, 
MRA-568, deposited by AF Cowman; DIO ACP (transit) -GFP, 
MRA-569, deposited by AF Cowman; W2, MRA-157, depos- 
ited by DE Kyle; 3D7HT-GFP, MRA-1029, deposited by A 
Talman and R Sinden; and NF54, MRA-1000, deposited by M 
Dowler, Walter Reed Army Institute of Research. P. falciparum 
strains were cultured in human O* erythrocytes as previously 
described.''' The DIO ACP-GFP transfectant strain was supple- 
mented with 100 nM pyrimethamine (Sigma-Aldrich, P7771).'^ 
Drug-treated parasites were cultured in the presence of final 
concentrations of 25 |JlM artemisinin (Sigma-Aldrich, 361593) 
and 100 nM chloroquine (Sigma-Aldrich, C6628). Parasites 
were visualized by fixing thin smears of cultures in methanol for 



1 min at room temperature then staining with Modified Giemsa 
Stain (Sigma-Aldrich, GS128). 

HMM search of autophagy-related Pfam domains 
Functional Pfam domains found in ATG gene-encoded pro- 
teins were examined in 11 genomes: P. falciparum (release ver- 
sion 8.1), P. vivax (release version S.l), P. knowlesi (release version 
S.l), P. yoelii (release version 8.1), P. chabaudi (release version 
8.1), P. berghei (release version 8.0), T. gondii (release version 7.1), 
Cryptosporidium hominis (release version 4.5), C. parvum (release 
version 4.5), Saccharomyces cerevisiae, and Homo sapiens. The 
Pfam HMM profiles were downloaded at http://pfam.janelia.org 
from their Pfam_ls HMM library version 26.0. Additional infor- 
mation regarding Pfam accession numbers are given in Table 55. 
HMM searches were performed, as previously described, in a 
series of threshold E values, increasing stringency from E value 
< 0.5 to E value < 0.05 (data not shown)." An E value of 0.05 
gave the most consistent results when compared with previously 
published results. 

Generation of PfATGS antibody 

PfATGS was amplified from cDNA using primers for- 
ward: CA ATCTAGAC CCGGGA ATAA TTTTGTTTAA 
CTTTAAGAAG and reverse: CAAGGATCCC 
CATGGGGGTA TAGGGGACAT ACCGCTGC, bearing 
the restriction enzyme sites Xbal/Xmal and BamHI/NcoI, 
respectively (underlined). The fragment was inserted into the 
pGS21a plasmid vector (GenScript, SD0121) using the Xbal 
and BamHI restriction sites. PfATCS-His^^ recombinant pro- 
tein was expressed in arctic bacterial cells, pulled down using 
Ni-NTA agarose beads (Qiagen, 30210), and then purified and 
concentrated by 10- and 30-kDa size exclusion centrifugation 
(Millipore, UFC901008, UFC90300S). Anti-PfATGS antisera 
were raised in rabbits and were then purified using protein A/ G 
(Thermo Fisher Scientific, Pierce Custom Antibody Service). 

Cloning of P. falciparum GFP-fused ATG8 

The DNA corresponding to the P. falciparum ATG8 ortho- 
log (PF3D7_1019900; former ID: PF10_0193) was obtained 
by PGR from cDNA using primers forward: TCCC CCCGGG 
ATGCCATCGC TTAAAGACG and reverse: ACGC GTCGAC 
CTGCAGTTAT CCTAGACAAC TCTCACAAC, bearing the 
restriction enzyme sites Xmal and Sall/PstI, respectively (under- 
lined). The fragment was cloned into the pCC-1 vector using 
the Xmal and Sail restriction sites. GFP was then amplified 
from the pARL/GFP vector (Courtesy of Dr Jude Przyboski, 
Philipps-Universitat, Marburg, Germany) using the primers for- 
ward: CCTG CTCGAG ATGAGTAAAG GAGAAGAACT and 
reverse: TCCC CCCGGG TTGTATAGTT CATCCATGCC, 
bearing the restriction enzyme sites Xhol and Xmal, respectively 
(underlined). The fragment was cloned into the N terminus 
of the ^CC-\tPfATG8 vector to make ^CC-\ I GVP-PfATGS. 
Using the primers forward: CAGG CTCGAG ACTAGTA TGA 
GTAAAGGAGA AGAACT and reverse: GGTC ACCGGT 
CCCGGGGTCG AC TTATCCTA GACAACTCTC AC, bear- 
ing the restriction enzyme sites Xhol/Spel and Agel/Xmal/Sall, 
respectively (underlined), the GFP-PfATGS fragment was ampli- 
fied and inserted into the pARL vector using Xhol and Agel 
restriction sites. The plasmid 'pAKLlGW-PfATGS was then used 
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for transient expression of N-terminally GFP-tagged ATG8 in 
P. falciparum. All restriction enzymes used were from New 
England Biolabs. 

Generation of GFP-PfATG8 transient expressing P. falci- 
parum cell lines 

Parasite transfection was performed as previously described. 
3D7 parasite cultures were synchronized with 5% sorbitol and 
ring stage parasites were used for electroporation with 40 to 50 |JLg 
of plasmid DNA. Each electroporation was performed in a 0.2 cm 
cuvette (Bio-Rad, 165-2086), with time constants ranging from 
11-14 ms. The following day, mature parasites were separated 
by magnetic purification and uninfected RBCs were electropor- 
ated with 40 to 50 |xg of plasmid DNA. All electroporations were 
pooled and 50 |xl of fresh RBC were added to the culture. The 
transfected cultures were supplemented with 5 nM WR99210 
(Courtesy of Jacobus Pharmaceuticals) for drug selection."^ 

Subcellular fractionation 

Parasite subcellular fractions were generated as previously 
described.''* Parasites were extracted from red blood cells using 
0.15% saponin (Acros, 41923-1000) in PBS (Thermo Fisher 
Scientific, BP399-1 [11.9 mM phosphate, 137 mM NaCl, 
2.7 mM KCl]) on ice for 10 min, followed by 3 ice-cold PBS 
washes. Parasite pellet fractions were incubated with cytoplasmic 
lysis buffer (0.65% Igepal [Sigma-Aldrich, 18896], 10 mM TRIS- 
HCl, pH 7.5, 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 2 mM 
AEBSF [Thermo Fisher Scientific, BP63550], Ix Roche protease 
inhibitor tablet [04693159001]) for 5 min followed by centrifu- 
gation at 2,000 rpm for 10 min. The supernatant fraction con- 
taining cytoplasmic proteins was removed, samples were washed 
with PBS, and then subsequently incubated with nuclear lysis 
buffer (10 mM HEPES, pH 79, 0.1 mM EGTA, 0.1 mM EDTA, 
1.5 mM MgCl,, 420 mM NaCl, 0.5 mM DTT, 25% glycerol, 
2 mM PMSF, Ix Roche protease inhibitor tablet) for 20 min, 
followed by centrifugation at 6,000 rpm for 10 min. The super- 
natant fraction containing nuclear proteins was removed, samples 
were washed with PBS, and then incubated with solubilization 
lysis buffer (2% Triton X-100, 50 mM TRIS-HCl, pH 7.5, 1 mM 
EDTA, 1 mM EGTA, 2 mM PMSF, Ix Roche protease inhibi- 
tor tablet) for 30 min and were then pelleted at 13,000 rpm for 
10 min.^' The supernatant fraction containing membrane proteins 
was collected. Protein concentration of samples was determined 
by the Bradford method.^" Equal amounts of protein from each 
fraction were loaded into a 16.5% Tris-Tricine gel, transferred to 
a PVDF membrane (Bio-Rad, 162-0177), and the membrane was 
stained with Ponceau S as a loading control before western blot 
analysis.^" Verification of fractionation was done with the DIO- 
ACP (transit) cell line, which expresses GFP in the cytoplasm, and 
fractions were immunoblotted with antibodies to GFP (Abeam, 
ab290) and histone IT3 (Abeam, abl791) to confirm their pres- 
ence in the cytosolic and nuclear fractions, respectively. 

Western blot analysis 

Parasite cultures were incubated in lysis buffer (1% Triton 
X-100, 50 mM TRIS-HCl, pH 7.5, 1 mM EDTA, 1 mM 
EGTA, 2 mM PMSF, Ix Roche protease inhibitor tablet) for 
30 min on ice and then mechanically lysed by needle and syringe. 
Samples were then centrifuged at 10,000 rpm for 15 min and 



the supernatant fraction was removed. Equal amounts of protein, 
10 |JLg, were used for immunoblotting. Regular protein separa- 
tions were performed using 16.5% Tris-Tricine SDS-PAGE gels 
(Bio-Rad, 456-3064). Primary antibodies anti-PfATG8 (Thermo 
Fisher Scientific, Pierce Custom Antibody Service) and anti-GFP 
were diluted 1:1,500 and 1:3,500, respectively. Secondary anti- 
body goat anti-rabbit IgG HRP conjugate (Bio-Rad, 170-5046) 
was incubated at a dilution of 1:25,000. 
Gametocyte induction 

Gametocytes were induced in the NF54, DIO ACP-GFP, 
and GFP-PfATG8 cell lines following a protocol adapted from 
Ifediba and Vanderberg.''' Parasites were synchronized by 5% 
sorbitol lysis and diluted the following day into 75-cm^ flasks 
to reach 0.5% parasitemia at a hematocrit of 8.3% (total volume 
of 15 ml). Cultures were stressed for 3 d by daily replacement 
of 10 ml of culture medium. Cultures with 5-10% parasitemia 
were then induced by increasing the medium to a final volume 
of 25 ml per flask. For the next 5 d, cultures were maintained by 
removing 10 ml of medium and adding 10 ml of fresh medium 
supplemented with 50 mM N-acetyl glucosamine (Alfa Aesar, 
AAAA13047) to extinguish asexual parasites, followed by feed- 
ing with regular medium. 

Fluorescence microscopy 

P. falciparum parasites were permeabilized with 0.015% sapo- 
nin, fixed with 4% paraformaldehyde, and blocked with IFA buf- 
fer (2% BSA, 0.05% Tween-20, 100 mM glycine, 3 mM EDTA, 
150 mM NaCl, Ix PBS). Cells were incubated with anti-PfATG8 
antibody (Thermo Fisher Scientific, Pierce Custom Antibody 
Service; 1:400) followed by anti-rabbit DyLight 550 (abeam 
ab98489; 1:500). DAPI (Invitrogen, D21490) was included in a 
wash at a final concentration of 5 jJLg/ml. Images were acquired 
using the Olympus BX40 epifluorescence microscope or the 
Leica SP5 confocal microscope. 

Generation of apicoplast-minus P. falciparum 

Following a similar protocol described by Yeh and DeRisi,^' 
parasite cultures were synchronized using 5% sorbitol and ring- 
stage parasites at 3% parasitemia were grown in 6-well plates 
(Costar, 3506). The antibiotic chloramphenicol (Sigma-Aldrich, 
C1919) was added to a final concentration of 100 |jlM followed 
by incubation for 48 h. After parasites reinvaded red blood cells, 
IPP (Isoprenoids LC, IPPOOlLi) was added to a final concentra- 
tion of 200 |jlM to rescue apicoplast-minus parasites. Parasites 
were grown for an additional generation before sampling cultures 
for immunofluorescence assays. 

P. falciparum starvation 

P. falciparum cultures were synchronized by 5% sorbitol. Early 
and late stage parasites were grown for 6, 24, and 72 h in serum- 
free complete medium or in starvation medium (0.043 mg/ml 
gentamicin [GIBCO, 15710-064], 38.5 mM HEPES, 0.18% 
sodium bicarbonate, 0.003 mM NaOH in RPMI 1640 without 
L-glutamine and Phenol Red [Corning, 17-105]). Regular media 
changes were performed every 24 h to remove metabolic waste. 

SYBR Green I-based fluorescence assay 

Bafilomycin A^ (LC Laboratories, B-1080) was diluted to 
10 fJiM and serially diluted by thirds into clear bottom 96-well 
plates (Costar, 3904). P. falciparum cultures were added at a 
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2.5% hematocrit with a 1% parasitemia to 96-welI plates and 
incubated for 72 h at 37 °C. Plates were frozen at -80 °C over- 
night. Lysis buffer (20 mM TRIS-HCl, 5 mM EDTA, 0.008% 
saponin, 0.08% Triton X-100, and 0.2 jjil/ml SYBR Green I dye 
(Invitrogen, S7585) was added to thawed plates and incubated 
at 37 °C. The plates were read using the Spectra Max Gemini 
EM reader (Molecular Devices). Data were then analyzed with 
SoftMax Pro v5 (Molecular Devices Software, Inc.).''^ 
Immunoprecipitation 

The parasite cell lines 3D7, GFP-PfATG8, and HT-GFP were 
synchronized with 5% sorbitol. To enrich for ATG8 vesicles, 
samples were incubated with 75 nM of bafilomycin for 3 h 
before harvest. Parasites with or without bafilomycin A^ treat- 
ment were extracted from red blood cells using 0.15% sapo- 
nin. Approximately 2 x 10'" cells (1 ml of parasite pellets) were 
homogenized by needle and syringe in a total volume of 3.5 ml 
of lysis buffer [50 mM TRIS-HCl, pH 7.5, 150 mM NaCl, 5 
mM EDTA, 1 mM PMSF, 1% Triton X-100, 5 M-M E-64 (AG 
Scientific, Inc., E-2030), Roche complete EDTA-free protease 
inhibitor cocktail] . Lysates were centrifuged at 13,000 rpm for 10 
min and the supernatant fraction was transferred to a clean tube. 
The supernatant fraction was precleared with PureProteome 
protein G magnetic beads (Millipore, LSKMAGG02) for 1.5 h. 
Precleared beads were removed and either anti-PfATG8 or anti- 
GFP (Abeam, ab290) antibodies were incubated for 6 h at 4 |xg/ 
ml before fresh beads were added for 30 min. Our custom anti- 
PfATG8 antibody was found to not efficiently precipitate endog- 
enous PfATG8 (samples immunoprecipitated by the anti-ATG8 
antibody were analyzed by MudPIT; however, the endogenous 
ATG8 bait [PF3D7_1019900] was either detected with, at best, a 
couple of peptides and less than 5 spectral counts, or not detected 
at all [data not shown]). Resin-containing immune complexes 
were washed twice for each buffer with wash buffers A (1% NP40 
[BDH Laboratory Supplies, 56009-2L], 1 mM EDTA, 1% BSA 
[Invitrogen, 15561-020] in PBS), B (buffer A plus 0.5 M NaCl), 
C (1% NP40, 1 mM EDTA in PBS), and D (1 mM EDTA in 
PBS) 2 min, followed by magnetic bead capture. Proteins were 
first eluted with 0.1 M glycine, pH 2.8 and neutralized with 100 
mM TRIS-HCl, pH 8.0, followed by a second elution using the 
same solution, a third elution using 10% 1,4-dioxane (Thermo 
Fisher Scientific, Dl 11-500), and a fourth elution with elution 
buffer (50 mM TRIS-HCl, pH 6.8, 100 mM DTT, 2% SDS) 
for 40 min at 65 °C. Eluates were pooled and precipitated with 
20% trichloroacetic acid (TCA). The resulting pellet fraction 
was washed once with 10% TCA and 4 times with cold acetone. 
As a control for the GFP-immunoprecipitation, 2 independent 
batches of unsynchronized 3D7 cells expressing HT-GFP were 
prepared in parallel to synchronized 3D7 cells expressing GFP- 
ATG8 and subjected to the same GFP-affinity purification 
protocol. As a control to overexpressing GFP-ATG8 in cells, 2 
independent replicates of 3D7 unsynchronized cells expressing 
GFP-ATG8 were also subjected to immunoprecipitation using a 
nonspecific antibody. 

MudPIT analysis of immunoprecipitated proteins 
TCA-precipitated proteins were solubilized in 100 mM TRIS- 
HCl, pH 8.5, 8 M Urea, 5 mM Tris(2-carboxylethyl)-phosphine 



hydrochloride (TRIS-HCl, Sigma-Aldrich, T3253; urea, Sigma- 
Aldrich, U1250; TCEP-HCl, Pierce, 20490) before chloroacet- 
amide (Sigma-Aldrich, C0267) was added to a final concentration 
of 10 mM. Endoproteinase Lys-C (Roche Applied Science, 
11047825001) was used at 1:100 w/w to digest proteins overnight 
at 37 °C. Samples were then brought to a final concentration of 
2 M urea with 100 mM TRIS-HCl, pH 8.5 and 2 mM CaCl, 
(VWR Scientific Products Inc., EM3000) before adding trypsin 
(VWR Scientific Products Inc., EM3000) (Promega, V5111) at 
1:100 w/w for a second overnight digestion at 37 °C. Formic acid 
(VWR Scientific Products Inc., JT0129-1) was added to a final 
concentration of 5% before samples were loaded onto a split-tri- 
ple-phase fused-silica microcapillary column (Deactivated fused 
silica tubing, 0.250 x 0.360 mm: Agilent Technologies, 160- 
2255-30; Deactivated fused silica tubing, 0.100 x 0.360 mm: 
Molex Connector Corp, 2000023),''' and placed in-line with a 
quaternary agilent 1100/1200 series HPLC (ThermoScientific) 
LTQ ion trap mass spectrometer equipped with a custom-made 
nanospray ionization source. Full MS spectra were recorded on 
the peptides over a 400 to 1,600 m/z range, followed by 5 tandem 
mass (MS/MS) events sequentially generated in a data-dependent 
manner on the first to fifth most intense ions selected from the 
full MS spectrum (at 35% collision energy) .^^ The dynamic 
exclusion was enabled for 120 s. A total of 3 technical repli- 
cates were acquired for each of the 4 types of samples analyzed 
(early trophozoite with and without bafilomycin Aj treatment 
and late schizont with or without bafilomycin A^ treatment). 
GFP-immunoprecipitations from 2 batches of asynchronized 
3D7-HT-GFP cells and immunoprecipitations in the absence of 
an antibody from 2 batches of asynchronized GFP-PfATG8 3D7 
cells were also analyzed as negative controls. 
Proteomics data analysis 

The MS/MS data sets were searched using SEQUEST against 
a protein database consisting of 5,538 P. falciparum (PlasmoDB 
rerelease 9.1) and 34,521 H. sapiens (NCBI 2012-08-27 release) 
proteins, complemented with 177 sequences from usual contami- 
nants (human keratins, IgGs, proteolytic enzymes). To estimate 
false-positive discovery rates (FDR), each sequence was random- 
ized (keeping amino acid composition and length the same) and 
the resulting "shuffled" sequences were added to the "normal" 
database (doubling its size) and searched at the same time. The 
GFP-PfATG8 sequence was added to this database. To limit the 
number of false-positive "shuffled" peptides passing the selection 
criteria, the filtering criteria were strict, notably enforcing that 
both peptide ends matched the enzyme used for the digestion. 
The average protein FDR was 2.6%, peptide FDR was 1.06%, 
and the spectral FDR was 0.40%.''' 

To address the question of whether any differences were 
observable between the samples and the controls, we used mul- 
tiple criteria to filter proteins associating with ATG8. First, 
proteins had to be detected reproducibly in at least one of the 
4 types of samples analyzed, i.e., present in at least 2 of the 3 rep- 
licates acquired per sample type. Next, proteins were either not 
detected in controls or with a 5-fold increase based on normal- 
ized spectral counts. We narrowed down the list of significant 
proteins further by requiring that the proteins be at least present 
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in 4 out of the 12 analyses. Finally, for proteins passing these cri- 
teria but detected in the controls, we calculated the fold changes 
between the respective purifications using a z-statistics-based 
statistical framework (update QSpec) for the significance analy- 
sis of differential expression.'^ The z-score was calculated using 
the web submission at http://www.nesvilab.org/qspec.php/. The 
spectral counts and the length of each protein were used as input 
for the QSpec software. We only kept proteins with the addi- 
tional z-score filter of > 1.7. 
Gene ontology analysis 

Enrichment of gene ontology (GO) terms within each sam- 
ple was analyzed using the software package "topGO" (writ- 
ten in R and maintained by the BioConductor project) .^^ For 
each GO domain (i.e., cellular component, biological process 
or molecular function), we compared the proteins identified by 
MudPIT to the full proteome of P. falciparum using the "clas- 
sic" algorithm in combination with a Fisher's exact test. In addi- 
tion, GO enrichment was analyzed among all ATG8-associated 
proteins identified by MudPIT, irrespective of time point or 



drug treatment. Gene ontology structures were visualized for 
GO terms with a P value < 0.01 using the R software pack- 
age "Rgraphviz".''^ To illustrate overlaps in GO terms between 
the different samples, a Venn diagram was drawn including all 
GO terms with a P value < 0.05 using the R software package 
"venneuler".'''* 
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